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ABSTRACT 



Qi^ I Aims. We aim to study the dust ejected by main-belt comet (MBC) (300163) 2006 VW139 to obtain information on the ejection 

, mechanism and the spectral properties of the object, to see if they are compatible with those of "normal" comets. 

• ■ Methods. Broad-band images in the g and r band as well as a low-resolution spectrum in the 0.35 - 0.9 fim region were obtained 

i-^ ■ with the GTC telescope (Roque de los Muchachos Observatory, La Palma, Spain). Images were analyzed to produce a color map 

and derive a lower limit of the absolute magnitude H. A Monte Carlo scattering model was used to derive dust properties such as 
mass loss rates and ejection velocities as a function of time. The reflectance spectrum was compared to that of the well-studied MBC 

■ 133P/Elst-Pizarro. The spectrum was also used to search for CN emission and to determine the upper limit of the CN production rate. 
^ ■ Results. The reflectance spectrum of 2006 VW139 is typical of a C-class asteroid, with a spectral slope 5'^,=0.5±1.0%/1000A. It is 

, similar to the spectrum of 133P and other MBCs. No CN emission is detected in the spectrum. A CN production rate upper limit of 

3.76 X 10^^ J"' is derived. Images show that the MBC present a narrow almost linear tail that extends up to 40.000 km in the anti-solar 
direction and more than 80.000 km in the direction of the object's orbital plane. The color of the tail is slightly redder than the Sun 
(S' between 3-6 %/lOOOA). The Monte Carlo dust tail model derived the mass loss rates and ejection velocity as a function of time, 
^ , and the results show that the activity onset occurs shortly after perihelion, and lasts about 100 days; the total ejected mass is ~ 2x10* 

CN. kg. 

■ Conclusions. The spectrum of 2006 VW139 suggests that it is not a "normal" comet. The spectrum is typical of the other observed 
' MBCs. Even if no CN emission is detected, the more likely activation mechanism is water-ice sublimation. Like other well studied 

MBCs, (300163) 2006 VW139 is likely a primitive C-class asteroid that has a water-ice subsurface depth reservoir that has recently 

■ been exposed to sunlight or to temperatures that produce enough heat to sublime the ice. 

Key words, minor planets, asteroids: individual: (300163) 2006 VW139 - comets: general - techniques: photometric - techniques: 
spectroscopic 
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1. Introduction and established an upper limit for the CN production rate of 

5^- ecw< 1-3x1024^-1. 

C3 ; Asteroid (300163) 2006 VW139 (hereafter VW139) is the eighth Dynamical simulations show that MBCs are extremely un- 

main-belt asteroid observed to be "ac tive", i.e., with a cometary- nkely to originate in the place where comets came from: the 

like dust tail or coma (Hsieh et al. |20irK With a semi-major trans- Neptu nian belt (TNB) and the Oort Cloud (e.g., Fernandez 

axisof a = 3.0516 AU, an eccentricity e=0.2011, and incli- gt al. l2UU2l l. They li kely f ormed in situ. This is the case for 

nation /=3.2392 °, VW139 fits into the category of main-belt VW139 (Hsieh et al. IM3, which could also be a member of 

comets (MBCs) fol lowin g the pure phenomenological definition an asteroid collisional family (Novakovic et al. 2012 ). A recent 

by Hsieh & Jewitt (|2006t. work suggests that some icy ti-ans-Neptunian objects (TNOs) 

. J j^-ir^-imr, 1 i j ■ might have been delivered to the asteroid belt during the Late 

Images and spectroscopy of VW 139 have been presented m ,, n 1 j ,j ■ , ^r,r,r, 1 . 

. , ° , Xi T,, , J 1 1 1- Heavy Bombardment (Levison et al. 2009 ), but even those sim- 

Hsieh et al. (120121 ). The authors reported that the obiect presents w r 1 , j 1 • i- ^- 1 . • i.-. 

. -1 1- J — . 1-11 r , 1 • J ulations fail to produce low-mclination, low-eccentricity orbits 

a tail aligned with the orbital plane of about 1 arcmin, and ^^^j^ those of^BCs 

a short anti-tail of ~ 10". They claimed that this implies that 

the observed dust was produced by the action of a long-lived. The asteroid 133P/Elst-Pizarro is the first discovered MBC 

sublimation-driven emission mechanism. They also reported that (Elst et al. 119961 ) and so far the best characterized. It has been 

there is no evidence of any CN emission in the object spectrum observed to be active at every perihelion passage around the or- 
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bital quadrant following perihelion, and to be inactive in the two 
quadrants around aphelion (Boehnardt et al. '1997'; Toth 120061 
Hsieh et al. 2006, 2009 1010, Bagnulo et al. 2010). This recur- 
ring activity supports the hypothesis that 133P has an ice reser- 
voir, which periodically sublimates and elevates dusty material 
into the coma and tail region. The existence of ice in MBCs is 
surprising, as they are close to the Sun, but recent thermal mod- 
els suggest that water ice can survive at greater subsurface depth, 
even in the asteroid belt (Schorghofer, 2008). 

The MBCs 133P and 176P/LINEAR are members of the 
Themis collisional family. These two MBCs share surface prop- 
erties with other Themis family asteroids (see Licandro et al. 
120071 1201 Ij which indicates that they have formed in situ. Water 
ice has recently been discovered on the surface of (24) Themis 
(Campins et al. lIOTOl Rivkin & Emery l20T0l ) and (65) Cybele 
(Licandro et al.. l201 II ). Water ice on the surface of (24) Themis 
likely comes from the interior, which supports the hypothesis 
that there is a water ice reservoir below the surface of some 
MBCs. The spectral properties of 133? and 176P are signifi- 
cantly different from those of "normal" comets, i.e., comets scat- 
tered from the Oort Cloud or the TNB, which argues against a 
cometary origin (Licandro et al. l201 i\ . 

Hsieh et al. (2012) also presented numerical simulations 
showing that the asteroid is dynamically stable for > 100 Myr 
and searched for a potential asteroid family around the object. 
"A cluster of 24 asteroids within a cutoff distance of 68 . At 
70 ms"\ this cluster merges with the Themis family, suggesting 
that it could be similar to the Beagle family". 

Understanding the origin of MBCs is crucial. If they are 
formed in situ and if their activity is caused by water ice sub- 
limation, there should be water ice in many asteroids and they 
must have formed in the main-belt region, presumably at a time 
when the snow line was at ~3 AU, during the first Myr of the 
solar system. If they are captured TNBs or Oort Cloud comets, 
the mechanisms that drove them to their present orbits need to 
be understood. The existence of ice in main-belt objects is sur- 
prising given their proximity to the Sun, and presents intriguing 
opportunities for constraining the temperature, composition, and 
structure of primitive asteroids and our protoplanetary disk. It is 
therefore critical for understanding the physical conditions and 
the mechanisms of planetary formation, and it also addresses the 
question of the origin of Earth's water: if the outer main belt has 
a large population of asteroids with ice, they could have con- 
tributed to the water on Earth. Finally, this indicates the extent 
and origin of volatiles in asteroids that could be used as resources 
for space exploration. 

Late in 2010 we started a target-of-opportunity program 
at the 10.4m Gran Telescopio Canarias (GTC), located at the 
El Roque de los Muchachos Observatory (La Palma, Canary 
Islands, Spain) to obtain images and, when possible, visible 
spectra of new discovered MBCs. VW139 is the second MBC 
observed within the frame of this program (see Moreno et al. 
12010 ). and in this paper we present the resulting data. Images 
and spectra together with details on data acquisition and reduc- 
tion are described in Section 2. An analysis of the data is pre- 
sented in Section 3. In Section 4 a Monte Carlo scattering dust 
model is applied to fit the observed dust tail and derive the dust 
properties, and to obtain information on the activation mecha- 
nism. Finally, the conclusions are presented in Section 5. 

2. Observations and Data Reduction 

Broad-band images and low -resolution spectroscopy of VW139 
were made in service mode on November 29, 2011 using the 



Date (UT) 


r(AU) 


A(AU) 


a{°) S uiipaC) Motionpf^(°) 


29/11/2011 


2.52 


1.67 


14.0 65.2 247.5 



Table 1. Observing geometry, r and A are the heliocentric and 
geocentric distance, a is the phase angle, SunpA and MotioripA 
are the position angles of the extended Sun-target radius vector 
and the negative of the target's heliocentric velocity vector. 



Optical System for Imaging and Low Resolution Integrated 
Spectroscopy (OSIRIS) camera-spectrograph (Cepa et al., 20001 
Cepa. I20T0I ) at the 10.4m Gran Telescopio Canarias (GTC), lo- 
cated at the El Roque de los Muchachos Observatory (La Palma, 
Canary Islands, Spain). The OSIRIS instrument consists of a 
mosaic of two Marconi CCD detectors, each with 2048x4096 
pixels and a total unvignetted field of view of 7.8'x7.8', giving 
a plate scale of 0.127 "/pix. To increase the signal-to-noise ra- 
tio (S/N) for our observations, the data were 2x2 binned, corre- 
sponding with the standard operation mode of the instrument. 
The observing geometry is presented in Tabl^U 

2.1. Imaging 

CCD images were collected under photometric conditions, dark 
moon, and with an average seeing of 0.9". Four consecutive se- 
ries of 5x30 sees in the Sloan g and r filters were obtained be- 
tween 21:05 and 22:29 UT. Bias correction, flat-fielding, and 
bad-pixel-masking were performed using standard procedures. 
The images were finally aligned on the MBC. 

Images were flux-calibrated using Sloan photometric stan- 
dards. Each group of five "on-object" aligned images were fi- 
nally median-averaged to obtain four combined images in g band 
and another four images in r band. In the first two g, r series the 
asteroid was too close to a very bright star, therefore only the 
images of the last two series were used for the analysis. 

A constant sky level, measured in image regions out of the 
observed tail, was subtracted from the images, and they were 
finally flux-calibrated in A/, i.e., the albedo multiplied by the 
filling factor of the dust in the coma (A'Hearn et al. ll995i Tozzi 
et al. 120071) . using the following formula 

A/ = 5.34 X 10" (^)' C, X l0-0-4(z.-^'\ (1) 

where r/, is the heliocentric distance in AU, dx is the detector 
pixel size in arcsecs, C, is the pixel signal in e"/s, and Zp and Mj 
are the zero points and the solar magnitude in the used filter. An 
example of one of the r-band-calibrated images is shown in Fig. 

2.2. Spectroscopy 

The visible spectrum of VW139 was obtained on the same night 
immediately after the images. We used the R300B and R300R 
grisms, with dispersions of 2.48 and 3.87 A/pixel, respectively. 
A 5" slit width was used, oriented toward the parallactic angle 
to minimize slit losses due to atmospheric dispersion. An image 
of the object in the slit is shown in Fig|2l Three spectra with an 
exposure time of 300 sec each were obtained with each grism, 
covering the 0.35-0.70 yum and 0.50-1.0 fim spectral ranges at 
an airmass of X = 1.1. 

Images of the spectra were bias- and flat-field-coiTected 
with lamp flats. The two-dimensional spectra were wavelength- 
calibrated with Xe+Ne+HgAr lamps. After the wavelength cal- 
ibration, sky background was subtracted. We fitted a first-order 
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Fig. 2. Final reflectance spectrum of VW139. We also show an 
image of the asteroid inside the slit. 



Fig.l. Af calibrated SLOAN r-band image of VW139. North 
is up, ast is to the left. The antisolar direction is at PA=65.2°. 
The image covers a region of 75000x75000 km^. This is the last 
image of the series acquired in r, i.e., the one with the near-by 
bright star at the maximum distance (visible in the E direction). 
The color look-up table is logarithmic, with black (left side of 
the color bar) corresponding to Af-Q and yellow (right side of 
thebar)toA/=l x 10"^ 



3. Data analysis 

3. 1 . Analysis of the images 

In Fig. [T] the object clearly presents a narrow and almost linear 
tail that extends up to 40.000 km in the anti-solar direction (po- 
sition angle PA~ 63 °) and more than 80.000 km in the direction 
of the object's orbit plane (PA~ 250°). This image is typical 
of comets when observed with the observer very close to the 
comet's orbital plane and, as we show in Sect. 4, implies a dust 
ejection that lasted several months. 

We first studied the color of the dust tail by computing the 
spectral slope S ' from the r-band by the g-band images using 



constant sky level using a section across the spatial direction, 
and selecting two regions close enough to the MBC, but care- 
fully avoiding any contribution from the observed tail. This pro- 
cedure was repeated across the spectral direction every 50 pix- 
els. After background subtraction, a one-dimensional spectrum 
of the optocenter of the MBC was extracted using a very nar- 
row extraction aperture of only + 3 pixels (the FWHM measured 
for the central part of the MBC is ~ 6 pixels) to minimize the 
dust contribution. To coiTect for telluric absorption and to ob- 
tain the relative reflectance, a G star from the Landonlt ( |1992| l 
list was observed using the same spectral configuration at an air- 
mass (X - 1.2) similar to that of the object, and was used as a 
solar analogue star. Each individual spectrum of the object was 
then divided by the corresponding spectra of the solar analog and 
then normalized to unity at 0.55 fj.m. The three individual spectra 
were averaged, obtaining a mean spectrum for each of the two 
grisms. Those spectra were finally merged using the common 
wavelength interval 0.60-0.70 jum. The final reflectance spec- 
trum in the 0.38-0.9 fj.m region is shown in Fig.|2l 

The MBC spectrum was also flux-calibrated using observa- 
tions of the standard star G191 (observed the same night) and 
standard procedures. The flux-calibrated spectrum of VW139 
was used to derive an upper limit for the CN production rate 
(see Sect. 3.2). 



S ' = (Fr/Fg * Fgs IF,s - l)/(6.228 - 4.844), (2) 

where Fr and Fg are the flux-calibrated r- and g images, and 
Fgs and F^s are the flux of the Sun in the g and r band. The 
Fgs /Frs ratio is computed by using the solar color (g-r)s=QA4-. 
The central wavelength of the g and r filters are 4844 and 6228 
A. 

The color of the dust in the tail is slightly redder than the 
Sun, with 5' between 6% and 3%/1000 A. 

The absolute magnitude of VW139 is H-16.1 according to 
the Minor Planet Center (MPC), and was computed when the 
asteroid was not active. We performed aperture photometry of 
the MBC using an aperture radius of 5" to determine a lower 
limit for the asteroid magnitude (asteroid + ejected dust in the 
aperture) and a lower limit for H on both filters. 

We applied this only to the last two combined images of 
each filter (because a bright star was too close to the asteroid 
in the first two images). The measured magnitudes are r=19.07 
and 19.14, and g= 19.60 and 19.63, and thus the correspond- 
ing absolute magnitude lower limits are Hr= 15.16+0.04 and 
//g=15.67+0.02 (assuming the typical G value of 0.15 of the 
lAU H,G system). Of course this is a lower limit, because the 
dust of the tail contributes to the asteroid luminosity. The re- 
sults are compatible with the value of H given by the MPC. 
Considering that we assumed G = 0.15, the uncertainties on the 
derived Hr and Hg lower limits are about 0.1-0.2 mag. Our Hr 

determination is compatible with the Hr - 15.1 15.0 + 0.1 

reported by Hsieh et al. (I2012I I between Nov. 05 and Dec. 04, 
2011. 
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3.2. Analysis of the spectrum 

The obtained reflectance spectrum of the optocenter of VW139 
is very similar to that of a C-type asteroid, almost neutral in the 
0.5-0.9 //m region (with a specti-al slope5'j,=0.5+1.0%/1000 A), 
and a UV drop below 0.55 fim. 

Because the object was active, the obtained spectrum is "con- 
taminated" by scattered light from the dust tail. The photometry 
reported in the previous section and the absolute magnitude re- 
ported in the MPC show that in an aperture of 5 arcsec, the flux 
from the dust is about twice that reflected by the asteroid sur- 
face. Hsieh et al. (2012) reported a value of 2.5. For that reason, 
and considering that the PSF of the asteroid is much steeper than 
that of the tail, we used a very small extraction aperture of 0.8 
arcsec to minimize the dust contribution. In this small aperture 
it is too difficult to evaluate how significant the contribution of 
the dust is, but in any case it is less than for a larger aperture. 
The difference in color between this spectrum {S '^,=0.5) and the 
color measured in the tail in the previous section (5 '=0.5) sug- 
gests that a large part of the light in this aperture comes from the 
asteroid surface. 

However, the color of the dust close to the asteroid is slightly 
redder than the Sun, as shown in the previous section. The spec- 
tra of the dust cometary coma are usually featureless and reddish 
in the wavelength range of the presented VW139 spectrum, but 
the spectrum of the asteroid itself might be bluish, as is typical 
of a B-type asteroid. A spectrum of VW139 when it is inactive 
is desiderable to accurately determine its taxonomical classifica- 
tion. 

A bluish to almost neutral spectral slope is atypical of active 
comets and in cometary nuclei (Kolokolova et al. i2004i Licandro 
etal.. 120071) . 

Hsieh et al. ( 12012!) also published a spectrum of VW139 ob- 
tained a week before our observations. To compare ours with 
their results, both spectra are shown in Fig. [3] The spectrum of 
Hsieh et al. (120121) . with a lower S/N, is very similar to ours be- 
low 6500 A and slightly redder at longer wavelengths. Hsieh et 
al. ( [2013 reported a spectral slope S ^=7.2%/1000 A, computed 
using the 0.4-0.9 fim spectral range. If, as is usually done by 
convention, we only consider the 0.5-0.9 fim region (to avoid 
the UV drop-off), the spectral slope is 5';,=5.4+1.0%/1000 A, 
still slightly redder than ours. 

The similarities between both spectra in the blue region (be- 
low 5500 A) is a clear indication that the differences are not 
caused byo problems related to differential extinction (the Hsieh 
et al. spectrum was obtained without the orientation toward the 
parallactic angle) and/or by slit losses, because, in that case, the 
blue region of the spectrum is the most affected. One possibil- 
ity is that because the Hsieh et al. spectrum was obtained with 
the slit aligned with the tail, maybe the aperture used to extract 
their spectrum was wider than the one we used (1.5"), and con- 
sequently they might have included a larger amount of dust. 

Licandro et al. (201 1 ) analyzed the nature of MBCs based on 
their spectra. In Fig.[3]our spectrum of VW139 is also compared 
to the spectrum of 133P/Elst-Pizarro published by Licandro et al. 
The spectra are very similar, therefore the conclusions about the 
nature of 133P in Licandro et al. (201 1 ) can be reproduced here 
for VW139: the spectrum of VW139 resembles that of asteroids 
of the C-complex and is significantly different from those of the 
large majority of observed cometary nuclei, typically D- or P- 
types. 

We also searched for any signature of gas emission that is 
typically observed in cometary comae, of which that of CN 
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Fig. 3. Comparison of the reflectance spectra of VW139 pre- 
sented in this paper (blue) and that from Hsieh et al. (2012) 
(green). We have also included in the plot the reflectance spec- 
trum of 133P/Elst-Pizarro (red) from Licandro et al. (2011b). 
The two spectra of VW139 are shifted vertically by 0.5 for clar- 
ity. 

(B^S^-X^S^; Av = 0) at -3800 A is the most prominent. The 
detection of gas emission would be a strong test to determine if 
the observed tail and coma are produced by ice sublimation as 
in normal comets. 

In Fig. ID we show a zoom of our flux-calibrated spectrum of 
VW139 in the wavelength region of this CN band and compare it 
with the spectrum of active Jupiter family comet 8P/Tuttle (Lara, 
L., personal communication). None of the molecular bands, in 
particular the CN emission, which is clearly seen in the spec- 
trum of 8P centered at 3870 A , are detected in the spectrum of 
VW139. Hsieh et al. (I2012I I also failed to detect CN emission in 
their lower S/N spectrum obtained a week before ours. 

An upper limit to the CN production rate is achieved by 
studying the flux measure in the 3830-3905 A range in an aper- 
ture of 40 pixels, i.e., 10.16 arcsec, meaning 12300 km at the 
comet distance. We also measured the continuum that borders 
the band at 3770-3815 A and 3910-3970 A by approximating 
the continuum contribution by interpolating the left- and right- 
hand continua. Any remaining flux in the 3830-3905 A region 
represents a Icr upper limit to the CN emission. 

We converted the emission band flux into column density 
using the g-factor from Schleicher (11957b at the asteroid helio- 
centric velocity of 2.124 km/s and distance of 2.524 AU. The 
resulting g-factor is g=3.217 x 10"' erg s"' molecule"'. 

To compute the gas production rate, we assumed the Haser 
modeling (Haser, 1957) with the CN parent velocity Vp scaled 
with r/, (vp - 0.86r^'''' km/s), and Vd - I km/s for CN itself. 
The CN parent scale length is 13.000 km, whereas the CN one 
is 210.000 km (AHearn et al. [T995] l. With this set of parame- 
ters, the Haser modeling provides a CN 1 cr production rate upper 
limit of 3.76 X 10^^ s-\ i.e. 3cr of 1.13 x 10^"* s'K 

Our determination is on the same order of magnitude as the 
one reported by Hsieh et al. (120121 1. who obtained a value of 
1.3 X 10'-* .r'. 

However, this does not mean that there is no gas emis- 
sion and/or that a reasonable upper limit of the production rate 
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Fig. 4. Visible spectrum of VW139 (solid line) in the region of 
the CN and C3 emission bands (in ergs cm"^ s"' A"') compared 
to the spectrum of active comet 8P/Tutle (dashed line) obtained 
on Jan. 1, 2008 (Lara, L. personal communication). The spec- 
trum of VW139 is averaged over an aperture of 40 pixels, i.e., 
10.16 arcsec, meaning 12300 km at the comet distance, and is 
multiplied by 50 for a better comparison with the spectrum of 
8P/Tuttle. 



of H2O (Qh^o) can be derived from it. A limit from Qh^o is 
usually determined assuming the average ratio Qcn/Qoh and 
QohIQh^o of Jupiter-family comets (e.g. Hsieh et al. 120121 
Licandro et al. l201 11 1. But, as noticed by Licandro et al. (1201 11 1, 
this is probably meaningless, because MBCs likely formed 
closer to the Sun than "normal" comets. This implies a higher 
temperature, therefore one would expect a much lower fraction 
of species that are more volatile than water ice with respect to 
water ice, i.e., QcnIQoh is probably much lower for the MBCs 
than for "normal" comets. 



4. The dust model 

We analyzed the asteroid images by a direct Monte Carlo dust 
tail model, which is based on previous works on cometary 
dust tail analysis (e.g. Moreno '2009', Moreno et al. 1201 Oi l, 
and on the characterization of the dust environment of comet 
67P/Churyumov-Gerasimenko before Rosetta's arrival in 2014 
(the so-called Granada model, see FuUe et al. l2010l l. 

Briefly, the code computes the trajectory of a large number 
of particles ejected from a cometary nucleus surface, which are 
exposed to solar gravity and radiation pressure. The gravity of 
the asteroid itself is neglected, which is a good approximation 
for small-sized objects. 

In an ice-sublimation-driven scenario, the particles are ac- 
celerated by gas drag to their terminal velocities, which are the 
input velocities considered in the model. Particle ejection by any 
other mechanism (e.g., a collision) can also be simulated. Once 
ejected, the particles describe a Keplerian trajectory around the 
Sun, whose orbital elements are computed from the terminal ve- 
locity and the ratio of the force exerted by the solar radiation 
pressure and the solar gravity (see Fulle 1989). This parameter 
can be expressed as yS = CprQpr/(2pr), where Cp,.=1.19x lO"-' 
kg m"^, and p is the particle density, assumed to be p=1000 kg 
m"^. For particle radii r > 0.25 fim, the radiation pressure co- 
efficient is Qpr ~ 1 (Burns et al. I1979I I. To compare the model 
results with the observed isophotes, for CPU and memory rea- 



sons, we rebinned the original image by 2x2 pixels so that its 
spatial resolution becomes 310 km pixel '. 

For the observation date, we computed the trajectories of a 
large number of dust particles and calculated their positions on 
the (A^, M) plane. M is the projected Sun-comet radius vector, 
and is perpendicular to M in the opposite half-plane with re- 
spect to the nucleus velocity vector. Then, we computed their 
contribution to the tail brightness. 

For a first approximation of the input parameters to impose 
on the dust model, we first obtained a syndyne-synchrone di- 
agram. Figure |5] shows this map with seven synchrones from 
-1-108 to -54 days to perihelion, and three syndynes correspond- 
ing to particle radii of 12, 25, and 100 /im. The M- and A^-axis 
are scaled independently to avoid label crowding. Based on this 
map, and since the observation was made about 130 days post- 
perihelion, we infer that the onset of the asteroid activity should 
have occurred around perihelion, lasting close to the observation 
date, as a first approximation. 

The synthetic tail brightness depends on the ejection veloc- 
ity law assumed, the particle size distribution, the dust mass loss 
rate, as well as on the geometric albedo of the particles. We as- 
sumed isotropic ejection. To constrain the number of free pa- 
rameters, the size distribution was assumed to be independent of 
time and to be characterized by a power-law with index a. The 
maximum size of the ejected particles is set by the escape veloc- 
ity. If //=16.1, and the albedo is ~0.05, the asteroid radius 
would be R -1.8 km. If the bulk density is p=1000 kg m"^, the 
escape velocity would be about 0.3 m s"', assuming a spherical 
object. Based on the syndyne-synchrone locus, we set the lower 
limit of particle sizes at 1 fim. Smaller particles ejected some 30 
days before the observation would be beyond the image limits. 

Then, we ran the code for a large number of input dust loss 
masses and particle velocities as functions of the heliocentric 
distance. The velocity was assumed to vary with yS as v oc yg'^^. 
The best fit is shown in Figure |6] whereas the dust loss rates 
and particle velocities are shown in Figure Q Our best-fit model 
used a power-law size distribution function with power index 
of -3.1, which is in the range of the estimates for other comets 
and MBCs. For the escape velocity, we found Vgsc=0.2 m s"' as 
the best-fit. With this value for the escape velocity, we derived 
from the model that the largest ejected particles forming the anti- 
sunward branch are about 200 fim in radius. If the escape veloc- 
ity is assumed to be higher than 0.2 m s"', the maximum size of 
those particles would be smaller, and the model would be unable 
to fit that anti-sunward branch. This is clear from the syndyne- 
synchrone network of Fig. |6] which indicates that this feature 
can be fitted by particles of size ~100 /zm or larger. On the other 
hand, if Vesc were slower than ~0.2 m s"' , it would be possible to 
fit the dust feature, because particles in excess of 200 /jm could 
be ejected. 

The total dust mass ejected derived from the model is ~ 
2x10'' kg, about 2x10"^ % of the total mass of the asteroid. The 
activity onset is shortly after perihelion, and lasts about 100 days 
(see Fig.|7]i. This mass loss peaks some 60 days after perihelion 
at 0.5 kg s and is correlated with the ejection velocity. This 
is compatible with the early observations of Hsieh et al. ( 1201 11 1, 
which describe the object as having a non-stellar PSF on August 
30, about 45 days after perihelion. On the other hand, the veloc- 
ities are consistent with expectations of cometary activity at 2.5 
AU from the Sun, which supports the assumption that water ice 
sublimation is the mechanism that ejected the observed dust. 
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M{kin) 

Fig. 5. Syndyne-synchrone diagram overimposed on a VW139 
image. The density of the particles is assumed at p=1000 kg 
m"^. The synchrones are shown as solid white lines, and cor- 
respond, in clockwise order, to H-108, +81, -1-54, +27, 0, -27, and 
-54 days to perihelion, as labeled. Syndynes are shown as white 
dashed lines, and correspond to 0.0012, 0.0025, 0.01, and 0.05 
cm, as labeled. The M- and A^-axis are scaled independently to 
avoid label crowding. The bright feature on the left centered at 
approximate coordinates (-5000,7000) corresponds to a a bright 
field star that can be seen in Fig.[T]placed near the middle of the 
left ordinate axis. 




60 



20 40 60 




Fig. 6. Left: median-averaged image of VW139 obtained on 
November 29, 2011, referred to the (N, M) plane. Center: mod- 
eled image. Right: Comparison of observed isophotes (black) 
and modeled isophotes (red). The contour levels are 5x10 
10 and 2x10 in solar disk intensity units. Each pixel cor- 
responds to 310 km. 



5. Discussion and conclusions 

The principal aim of studying MBCs is, as explained in the in- 
troduction, to understand their activation mechanism and origin. 
Are they comets, or asteroids that, under certain circumstances, 
eject dust? 

Images of the ejected dust can help us to understand what the 
dust ejection mechanism is. The images of 2006 VW139 obtained 
on November 29, 201 1 with the OSIRIS camera-spectrograph at 
the 10.4m GTC (El Roque de los Muchachos Observatory, La 
Palma, Canary Islands, Spain) show an object that presents a 
narrow, almost linear tail that extends up to 40.000 km in the 
anti-solar direction and more than 80.000 km in the direction of 
the object's orbit plane. This tail is typical of comets observed 
very close to the comet's orbital plane. Moreover, the color map 
shows that the tail in the anti-solar direction is significantly red- 
der, which indicates that the larger particles concentrate in this 
region, as expected for an object that presents a comet-like ac- 
tivity. 

Using a direct Monte Carlo dust tail model (e.g. Moreno 
120091 Moreno et al. 120 1 Oi l, we fitted the images and derived the 
mass loss rates and ejection velocity as a function of time. We 
first concluded that the observed tail is not compatible with a 
short-lasting event like a collision or an outburst. The activity on- 
set occurs shortly after perihelion, and lasts about 100 days. All 
evidence favors a comet-like ejection mechanism, which means 
dust ejected by the sublimation of volatiles. 

We were also able to derive other properties of the ejected 
dust and the object itself. From the dust model, and to fit the 
sunward branch intensities, the best-fit escape velocity derived 
was vvsc=0.2 m s"'. 

On the other hand, the mass loss peaks some 60 days after 
perihelion at 0.5 kg s ', and is correlated with the ejection ve- 




Days to perihelion 

Fig. 7. Derived mass loss rate (solid line, left vertical axis) and 
ejection velocity (dashed line, right vertical axis), as a function 
of time. The ejection velocity refers to a 0.01 cm particle. The 
dotted line indicates the assumed escape velocity. 



locity. This is compatible with the early observations of Hsieh et 
al. (1201 11) . which describe the object as having a non-stellar PSF 
(indicative of a dust cloud around the object), on August 30 (45 
days after perihelion). The derived velocities are also consistent 
with expectations of cometary activity at 2.5 AU from the Sun. 
The total mass ejected is ~ 2x10*' kg. 

The spectrum of the object can be compared with spectra of 
comets and primitive asteroids to determine which group the ob- 
ject belongs to, and to search for gas species in the coma (which 
would confirm of a comet-like activation mechanism). The re- 
flectance spectrum of the optocenter of VW139 is typical of a 
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C-type asteroid, almost neutral in the 0.5-0.9 fwa region (with 
a spectral slope 5'5,=0.5±1.0%/1000A), and a UV drop below 
0.55 jum. A bluish to almost neutral slope in the spectrum is 
atypical of active comets but typical of MBCs (Licandro et al. 
1201 I I). Accordingly, VW139 presents a spectrum typical of the 
asteroids with similar semi-major axis and different from those 
of comet nuclei. This supports an in situ origin, as for dynamical 
simulations. 

We finally searched for the CN (B^'L^-X^'Z^; Av = 0) sig- 
nature of gas emission at ~3800A that is typically observed 
in cometary comae, finding no evidence of this feature in our 
spectrum. We derived an upper limit of the CN production rate 
QcN - 3.76 X 10^^ s Unfortunately, the direct determination 
of GhjO is very difficult, and to derive an upper limit of Qh^o 
based on the upper limit of Qcn, we had to assume values of the 
QcnIQoh and QohIQhiO ratios that are typical of comets. This 
last assumption might not be correct for MBCs, as explained in 
Sec. [321 

In conclusion, (300163) 2006 VWng, like other well-studied 
MBCs (see the case of 133P/Elst-Pizarro) is likely a primitive C- 
type asteroid that has a water-ice reservoir at subsurface depth. 
Owing to still uncertain mechanisms (probably a collision that 
eroded its surface), this reservoir has recently been exposed to 
sunlight or to temperatures that produce enough heat to sublime 
the ice and drag dust from the surface. 
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